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A Theory for Fluctuations in Stock Prices and Valuation of their Options 
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A new theory for pricing options of a stock is presented. It is based on the assumption that 
while successive variations in return are uncorrelated, the frequency with which a stock is traded 
depends on the value of the return. The solution to the Fokker-Planck equation is shown to be 
an asymmetric exponential distribution, similar to those observed in intra-day currency markets. 
The "volatility smile," used by traders to correct the Black-Scholes pricing is shown to provide an 
alternative mechanism to implement the new options pricing formulae derived from our theory. 
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Although options contracts have been in use as far back 
as the reign of Hammurabi in Babylon [1] , the first use- 
ful theoretical analysis for their valuation was presented 
relatively recently [2,3] . In the Black-Scholes theory, vari- 
ation in the price S{t) of a stock are investigated using 
the "return" x{t) — log[S{t)/ Sq] [4], where 5*0 is a "con- 
sensus" value of the stock at the time [t = 0) that the 
option is purchased [5] . 5*0 is typically set to be the price 
S{t = 0). The expected growth rate b for S{t) satisfying 



(Sit)) = Soe 



bt 



(1) 



may differ from the interest rate r on funds borrowed to 
purchase the option. It is further assumed that succes- 
sive random fluctuation in returns are independent and 
identically distributed, and hence (by central limit the- 
orem) that x{t) lies on a normal distribution for suffi- 
ciently large t [6]. A European call (i.e., an option to 
purchase a stock at a "strike" price K), is valued by 
its expected profit at expiration; i.e., CssiS, K,t) = 
exp{-rt) dS{S - K)fLNiS,t) [3]. In the Black- 
Scholes theory, the distribution of stock prices fLN{S,t) 
is log-normal and it can be shown that [3] 

CBs{K,So,t) = 5oe('-'^)*AA(d+) -Xe-'^*AA(d_). (2) 

Here J\f(x) denotes the cumulative normal distribution, 

(3) 



, bt + log{So/K) aVi 
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and a is referred to as the volatility of the return. 
Similarly, a European put (i.e., an option to sell a 
stock at a price K), is valued using PBs{S,K,t) — 
exp{~rt) dS{K — S)Jln{S, t), and is given by 

PBs{K,So,t) = Ke-^'W{~d^) - 5oe(''-'')*AA(-d+). (4) 

Recent investigations of bond and foreign exchange 
markets have clearly shown that the distribution of re- 
turns (for a fixed delay) deviates significantly from a nor- 
mal distribution, especially far from the mean [7-12]. In 
fact, to a very good approximation, intra-day currency 



fluctuations lie on an an asymmetric exponential dis- 
tribution, see Figure 1 [13]. It is also known that cur- 
rency traders do not assign values for options according 
to Eqns. (2) and (4). The most common correction is 
the use of a heuristically deduced "volatility smile"; i.e., 
a strike price dependent expression for a (which contra- 
dicts the basis of the Black-Scholes theory). In this Let- 
ter, we argue that the use of a volatility smile can account 
for changes in the value of options due to deviation of the 
distribution of returns from normality. This conclusion 
is reached via a new theory that is constructed on the 
assumption that although successive events are uncor- 
related, the step size of the corresponding random walk 
depends on the values of the return x{t) and time t. 
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FIG. 1. The relative frequency of intra-day returns (delay 
of 4 hours) for U.S. bonds over a period of 600 consecutive 
trading days beginning in March, 1985. For returns shown, it 
is approximately an asymmetric exponential distribution. 

Analyses of financial markets have provided evidence 
for the assertion that successive variations of the return 
are uncorrelated [11,14]. It implies that the distribution 
function W{x, t) for returns x{t) satisfies a Fokker-Planck 
equation [15,6] 

^\iDW), (5) 
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D = D{x,t) being the diffusion coefficient. The drift 
rate B is assumed to be constant for intra-day fluctua- 



tions [16] . Following conclusions from studies of financial 

markets [7,11], wo limit our considerations to solutions 
that have a scaling form, 



W{x,t) = ^F{u). 



(6) 



Here u = x/f, and t] is referred to as the "drift ex- 
ponent." Notice that the pre-factor t~'^ is introduced in 
order that W{x, t) is normalized for all t. 

Consider first the case B = 0. Financial markets ex- 
hibit the following behavior; a stock whose price S{t) 
deviates significantly (up or down) from 5*0 is traded at 
a higher frequency; i.e., the local diffusion rate D{x, t) is 
enhanced. This effect is further amplified if the variation 
occurs during a shorter time interval. These observations 
are quantified in an assumption that D(x, t) is a bi-linear 
function of u, 

D{x,t) = - ju)e{u) + + i^u)e{u). (7) 

Here Q{u) denotes the Heaviside 9-function and Q{u) = 
1 — 0(u). As will become apparent momentarily, for 
W{x, t) to be asymmetric, it is necessary for the parame- 
ters 7 and f to be different. Notice that for larger values 
of u, nonlinear corrections to D{x, t) may be required. 

These assumption on the absence of correlations be- 
tween successive movements of the return, and the scal- 
ing hypotheses (6) and (7) provide a unique value for 77. 
This can be seen from Eqn. (5), which simplifies to 

Consequently 77 = 1/2, in agreement with conclusions 
from the analysis presented in Ref. [7] which show that 
the standard deviation of W{x,t) for the S&P500 time 
series increases approximately as a{t) ~ t'^'^^. 

Next, substituting the diffusion coefficient (7) in 
Eqn. (8) for a; > gives, 



1 + vu 



F"{u)+(u+'^F'{u) + F{u) = Q. (9) 



Writing F(u) = E„=o«"« [17], it is found that 
F{u) ~ exp{—vu) is one of the solutions [18]. Combining 
it with the corresponding analysis for a; < provides the 
solution 



W{x,t) = ^eT"e(w) + ^e" 



^e(u) 



(10) 



to the Fokker-Planck equation. Two conditions are re- 
quired to evaluate the coefficients A and B. Normalizing 
W{x,t) gives 



7 v 



(11) 



The second condition imposed is the continuity of the 

probability current j{x,t) = —{DW)x, without which 
probability would accumulate at the origin. Using 
e^{u) = e(u), Quiu) = 6{u), e„(u) = -S{u) and 
f{u)S{u) = f{0)6{u), 
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Hence, continuity of j{x,t) at « = implies that 



A 



B 

~2 ■ 



(12) 



(13) 



Thus A = 72/(7 + u) and B = i''^/{-y + v). Notice 
that W{x^t) contains a discontinuity at the origin of 
(7 — v)l\ft. Notice also that the mean value (a;(t)) = 
for all t and that the variance is "21 /^v; in particular, 
W{x^t = 0) = 5{x)^ as required by the initial condition. 

We make two further observations. If D{x, t) is trans- 
lated by A, then the corresponding solution to the 
Fokker-Planck equation is W{x — A,t); however, the re- 
quirement {x{t = 0)) =0 forces A = 0. Second, the 
solution for the general case (i.e., B ^ 0) is obtained by 
Galilean transformations of D{x,t) and W{x,t). This 
can be achieved by replacing u by (x — Bt) / \/t. For this 
case {x{f)) = Bt. 

These conclusions were tested by investigating the fol- 
lowing two types of random walks (for the case B = 0). 
In the first, fixed-step-size walk (i.e., each step has a unit 
magnitude), the time interval for a step from a location 
x{t) is chosen to be 7^/(1 — ju) if m < and to be 
v^/ (1 -|- uu) if u > 0. In the second, fixed-step-time walk 
(i.e., each step takes a unit time) the size of a step is 
chosen to be ■\/(l — ju)/j or {1 + vu)/v depending on 
whether u < or u > 0. In each case the motion begins 
at the origin, and the direction of each step is chosen 
randomly with equal probability. The effective diffusion 
coefficient for these random walks is given by Eqn. (7) [6]. 

Figure 2(a) shows the histogram of positions for 5 mil- 
lion fixed-step-time walks of 256 steps, which is clearly 
consistent with W{x,t) given by Eqn. (10). The mean 
value of the power spectra for these random walks, shown 
in Figure 2(b), exhibits a fc"^ decay; a similar behav- 
ior has been reported in an analysis of financial mar- 
kets [7,19]. 

One final point needs to be clarified prior to deriving 
formulae for valuation of options. Condition (1) implies 
that 



(14) 



(-1 + ^^t){v-^^t) 



Since this is not an identity, the equality cannot be valid 
for all t; hence a growth rate given by (1) is not consistent 
with the new theory. We propose to replace it by 



{x{t)) 



(15) 



This last condition imposed on our theory imphes that 
the parameter B in the Fokker-Planck equation is b. 
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FIG. 2. (a) The relative frequency of final locations for 5 
million fixed-step-time walks of 256 steps. In this example 
7 = 10.0 and v = 15.0. The dashed lines show the distri- 
bution (10). (b) The mean {S{k)) of the power spectrum 
for these random walks exhibits a decay, similar to that 
observed in analyses of financial markets. 

Now we are in a position to re-derive the formulae for 
pricing of options. Consider first, a European call to be 
exercised at time t. Pricing the option by its expected 
profit, it is found that 



s''*G' - K 



KC 



1 + Vt 
l-Jt 



z < 1 

z > 1 



(16) 



where z = {K / SQ)exp{—bt) and ( = \/i/(7 4^ v). 

A European put can be valued using the definition 
P{SQ,K,t) = e-"-* dS{K - S)f{S,t). When the re- 
turn is distributed according to Eqn. (10), 

K - e^^G + KC- ^ 



Pe'' 



= KC- 



1 



1- 



Vt 



z > 1 

Z<1 



(17) 



Eqns. (16) and (17) can be used to justify the use 
of volatility smile curves to correct Black-Scholes pric- 



ing formulae. Suppose the "true" distribution of stock 
returns is (10), which decays slower than a Gaussian. 
Larger deviations of x{t) will occur more frequently than 
predicted by the Gaussian distribution under this sce- 
nario (and will be noted by an observant trader). As 
a result, these options will be valued higher than the 
pricing given by the Black-Scholes theory. For a given 
strike price, this enhancement can be accounted for by 
increasing the "effective" volatility (used in the Black- 
Scholes formulae). As an example, consider a distribu- 
tion W{x,t) with 7 = 15.0 and ly = 10.0. (The an- 
nualized volatility for this distribution is 11.5%.) The 
"correct" valuation for options is assumed to be given by 
Eqns. (16) and (17). We assign an effective volatility 
Vef / at a strike price K by equating the sum of prices of 
a call and a put estimated from the Black-Scholes theory 
(i.e., Eqns. (2) and (4)) with the corresponding sum for 
the new theory; i.e., Veff{K) is chosen so that 

PssiSo, K, t) -f Cbs{So. K, t) = P(5o, K, t) + B{So,K, t). 

(18) 



The function Veff{K) (for 5*0 = 100, t = 0.5 years, 
b = 10%, and r arbitrary), shown in Figure 3, resem- 
bles volatility smile curves used in financial trading. 




FIG. 3. The effective volatility evaluated using Eqn. (18) 
for several values of the strike price K, in the neighborhood 
of the underlying stock price of 100 units. 7 and 1/ are 10.0 
and 15.0 respectively, and the time to expiration of the op- 
tion is half a year. The dashed line shows a quartic fit to 
the data, and the dotted line the standard deviation for the 
distribution. 

We conclude by elaborating on a couple of issues raised 
by our work. The consensus value of a stock depends 
on many factors such as its historical performance and 
the market's expectations for its future prospects. In 
the Black-Scholes theory, it is used to define the returns. 
Typically, 5*0 is chosen to be the price of the stock at 
t = 0. If a different value is used for 5*0 (with a suitable 
modification of Eqn. (1)), the only effect is a uniform 
shift of W{x,t); the pricing formulae remain unchanged. 



(Another way to state this is that the Langevin equa- 
tion [6] is independent of 5*0.) In contrast, Sq plays a 
unique role in our theory. It is the value of the stock at 
which the diffusion coefficient reaches a minimum. Any 
deviation of the stock price from this value is reflected 
in an increase in the magnitude of its fluctuations. (The 
Langevin equation depends on through the diffusion 
coefficient D{x,t).) 
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FIG. 4. A solution F{u) of the Fokker-Planck equation 
with a diffusion coefficient that includes a quadratic term in 
u. For this example, "/ = v — 1.0, and e — —0.2. 

The second issue concerns possible higher order cor- 
rections to the bilinear form (7) for the diffusion coeffi- 
cient. For example, if a quadratic term 7~'^eu^ is added 
to D{x,t) for a; > 0. Then, the ordinary differential 
equation satisfied by F{u) changes to 

{l + iyu + €u^)F"{u) + {2v + Aeu + v'^u)F'{u) 

+ {2e + v^)F{u) = Q. (19) 

The solution for the symmetric case (7 = v) to this equa- 
tion that satisfies the initial condition F'{Q+) = —i'F{0) 
is shown in Figure 4. F{u) can be shown to decay as 
a power law for large u. Observations of such "fat- 
tails" have been made on distributions of inter-day re- 
turns [11,7]. In spite of such possible corrections to the 
distribution W{x, t), we believe that the theory presented 
here is still relevant in determining valuation of options. 
This is because, current frequency of trades (and not 
long time behavior of markets) is the information imme- 
diately available to traders, and on which their decisions 
are based. Hence, one may expect that distributions of 
intra-day variations in the return will determine the val- 
uation of options. 
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